Applications of Little's Law to stochastic models of gene expression 
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The intrinsic stochasticity of gene expression can lead to large variations in protein levels across a 
population of cells. To explain this variability, different sources of mRNA fluctuations ('Poisson' and 
'Telegraph' processes) have been proposed in stochastic models of gene expression. Both Poisson 
and Telegraph scenario models explain experimental observations of noise in protein levels in terms 
of 'bursts' of protein expression. Correspondingly, there is considerable interest in establishing 
relations between burst and steady-state protein distributions for general stochastic models of gene 
expression. In this work, we address this issue by considering a mapping between stochastic models 
of gene expression and problems of interest in queueing theory. By applying a general theorem 
from queueing theory. Little's Law, we derive exact relations which connect burst and steady-state 
distribution means for models with arbitrary waiting-time distributions for arrival and degradation 
of mRNAs and proteins. The derived relations have implications for approaches to quantify the 
degree of transcriptional bursting and hence to discriminate between different sources of intrinsic 
noise in gene expression. To illustrate this, we consider a model for regulation of protein expression 
bursts by small RNAs. For a broad range of parameters, we derive analytical expressions (validated 
by stochastic simulations) for the mean protein levels as the levels of regulatory small RNAs are 
varied. The results obtained show that the degree of transcriptional bursting can, in principle, be 
determined from changes in mean steady-state protein levels for general stochastic models of gene 
expression. 
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INTRODUCTION 

The intrinsic stochasticity of biochemical reactions in- 
volved in gene expression can give rise to large variations 
in protein levels across an isogenic population of cells 
[l|, y. Variations in protein levels, in turn, can give rise 
to phenotypic heterogeneity and non-genetic individual- 
ity in a population of cells [3|. The potential benefits 
of such phenotypic heterogeneity have been discussed for 
diverse systems [j| . Correspondingly there has been con- 
siderable effort focusing on uncovering molecular mecha- 
nisms which drive gene expression 'noise' as a source of 
phenotypic heterogeneity. 

High variablity in protein levels has generallybeen at- 
tributed to fluctuations in mRNA synthesis [5|-l7|. To 
elucidate the source of mRNA fluctuations, two distinct 
models have been proposed Q. In one case, mRNA syn- 
thesis is modeled as a Poisson process (Poisson scenario) 
and the high variability in protein levels is related to low 
abundance and infrequent synthesis of mRNAs [Qj. In 
the other case, fluctuations are primarily driven by the 
slow kinetics of promoter switching between active and 
inactive states (Telegraph scenario) with mRNA synthe- 
sis occuring only during the active stage [5|, llO| . Recent 
work has further generalized these models for gene ex- 
pression to include the effects of processes that can give 
rise to 'gestation' and 'senescence' periods for mRNA 
birth and decay ll|. These terms derive from the ob- 



servation that both creation and degradation of cellular 
macromolecules (mRNAs/proteins) often involve multi- 
ple biochemical steps. Correspondingly, the waiting-time 
distributions for these processes are more general than 
simple exponential distributions which are characteristic 



of single-step Poisson processes. Since the observed noise 
in protein levels can include contributions from different 
sources (e.g. transcriptional bursting as well as gestation 
and senescence) most single-cell measurements of steady- 
state protein distributions cannot be used to determine 
the source of fluctuations in mRNA synthesis. 

Recent studies have determined the variation of noise 
in protein expression as a function of mean protein abun- 
dance for several genes [6|, |7| . The observed scaling re- 
lationship is consistent with both the Poisson and Tele- 
graph scenario models in the limit that protein produc- 
tion occurs in infrequent random bursts [6|, [l^ . Further- 
more, advances in single-molecule techniques have led to 
studies monitoring real-time synthesis of proteins in sin- 
gle cells 13[ |l4 1 . Protein expression was indeed seen to 
occur in random bursts with a mean separation between 
bursts that is large compared to typical mRNA lifetimes 
13|,|l4|. Given these features, it is of interest to consider 
whether quantification of protein burst distributions can 
discriminate between the two scenarios for mRNA syn- 
thesis. For the Poisson scenario, the observed burst is a 
consequence of translation from a single mRNA, whereas 
for the Telegraph scenario, it is produced from a ran- 
dom burst of mRNAs synthesized when the promoter is 
in the active state. While the underlying mRNA burst 
distributions can thus be distinct for the two scenarios, 
it can be shown that observations of protein burst dis- 
tribution do not uniquely identify the underlying mRNA 
burst distribution |15| . 

Given that protein steady-state and burst distributions 
cannot discriminate between the Poisson and Telegraph 
scenarios, it has been argued that dynamic measure- 
ments of the number of mRNAs in single cells are needed. 



Such methods have indeed been developed in recent years 
[3, ll^: ll6| , and have been used to quantify the degree of 
transcriptional bursting. In this context, it would be of 
interest to derive equations relating burst and steady- 
state distribution means for both mRNAs and proteins. 
Such relations can provide useful checks for experimental 
approaches for measuring mRNA/protein burst distribu- 
tion means. Furthermore they can also suggest alterna- 
tive approaches which allow inference of the underlying 
mRNA burst distribution. This work focuses on deriv- 
ing such relations between the means of mRNA/protein 
burst and steady-state distributions and exploring their 
consequences for approaches to quantify the degree of 
transcriptional bursting. 

In this paper, we consider a mapping between general 
stochastic models of gene expression [ll| and problems of 
interest in queueing theory. By applying a general theo- 
rem from queueing theory, Little's Law, we derive exact 
relations connecting mRNA/protein burst and steady- 
state distribution means for stochastic models of gene 
expression with arbitrary waiting-time distributions for 
arrival and degradation of mRNAs and proteins. Fur- 
thermore the derived relations can be used to show how 
mRNA burst distributions can be inferred from measure- 
ments of mean protein levels by introducing an additional 
interaction in the reaction scheme. Specifically, we con- 
sider a reaction scheme that includes interaction between 
mRNAs and regulatory genes called small RNAs. In bac- 
teria, small RNAs have been studied extensively in recent 
years 17[ in part due to the critical roles they play in 



cellular post-transcriptional regulation in response to en- 
vironmental changes. The results derived in this work, 
besides the potential applications for quantifying the 
degree of transcriptional bursting, also provide insight 
into small-RNA based regulation for specific parameter 
ranges. 



MODEL AND RESULTS 
Connecting burst and steady-state means 

We begin by considering the minimal reaction scheme 
for translation from mRNAs 
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A single burst corresponds to proteins produced from the 
underlying mRNA burst distribution until decay of the 
last mRNA. For the Poisson process, mRNA transcrip- 
tion occurs with constant probability per unit time km- 
On the other hand, for the Telegraph process, mRNA 
transcription occurs with constant rate km only when the 
DNA is in the active(ON) state; once it transitions from 
the ON state to the inactive OFF state (with rate a) no 
mRNA transcription can occur until it transitions back 



from the OFF state to the ON state (with rate /3). Note 
that in the limit /3 << a, mRNAs will be produced in 
infrequent bursts. Thus the Poisson process gives rise to 
only protein bursts, whereas the Telegraph process gives 
rise to both mRNA and protein bursts. It is interest- 
ing to note that the mRNA burst distribution for both 
the Poisson and Telegraph scenarios can be represented 
by the conditional geometric distribution; specifically by 
considering bursts conditional on production of at least 
one mRNA [l5|. This can be understood as follows: the 
mRNA burst distribution is the number of mRNAs pro- 
duced in the active state before transition to the inactive 
state. Let us take the initial condition to correspond to 
the transition from the inactive (OFF) state to the active 
(ON) state i.e. at t = the DNA has just transitioned 
to the ON state. Now the next reaction that can occur 
either results in the production of a mRNA (with rate 
km) or a transtion to the OFF state (with rate a). The 
probablity of the next reaction being protein production 
whereas the probability that it is a transition to 
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the Off state is , ", . If we set p — , ", , the number 
of mRNAs produced {m) before the transtion to the OFF 
state, conditional on the production of atleast 1 mRNA, 
is given by 



TT^{m) = {l-p)"'-^p, m>l, 

^m(0)=0. 



(2) 
(3) 



For the conditional geometric distribution given above, 
the mean is given by 

(4) 
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For the Poisson scenario (p = 1), a single mRNA is 
produced per burst, which corresponds to the conditional 
geometric distribution with mean ?7^f, = 1. The Telegraph 
scenario also gives rise to a conditional geometric distri- 
bution for mRNA bursts, but with mean miy > 1. Thus, 
determination of the degree of transcriptional bursting 
(mti) can discriminate between the Poisson and Telegraph 
scenarios for intrinsic noise in gene expression. 

The general model for gene expression that we analyze 
is as follows. Bursts of protein expression result due to 
translation from the underlying mRNA burst, which has 
a conditional geometric distribution with mean m;,. The 
number of proteins produced from different mRNAs are 
taken to be independent random variables. The decay 
time for mRNAs and proteins is assumed to be drawn 
from arbitrary waiting-time distributions with means Tm 
and Tp respectively. Likewise, the waiting-time distri- 
bution between consecutive bursts is a random variable 
drawn from an arbitrary distribution with mean T{,. Cor- 
respondingly, the average arrival rate for bursts is given 
by fcf, = — . Since the waiting-time distributions are arbi- 
trary, effects due to gestation and sensecence of mRNAs 
and proteins 11| are included. For this setup, we will de- 
rive analytical relations which can be used to determine 



nib and thereby to quantify the degree of transcriptional 
bursting. 

We begin with the observation that the processes con- 
sidered in the above model have exact analogs in prob- 
lems of interest in queueing theory. For example, the 
creation of proteins corresponds to the arrival of cus- 
tomers in queueing models 18'. On the other hand, the 
service-time distribution corresponds to the waiting-time 
for the customer to depart the system, making it the ana- 
log of the waiting-time distribution for degradation of 
proteins. Given that degradation of each mRNA/protein 
is independent of other mRNAs/proteins in the system, 
the mapping corresponds to queueing systems with in- 
finite servers. This can be seen as follows. In infinite 
server queues, since the number of servers is unlimited, 
each customer is associated with a server immeduately 
upon arrival. This effectively implies that each customer 
is served independently of the others, which for the gene 
expression model is equivalent to the assumption that 
mRNAs/proteins are degraded independently. 

A general theorem from queueing theory, Little's Law 
|18| , states that the average number of customers in the 
system (L), the mean arrival rate (A) and the mean wait- 
ing time of a customer in the system {W) are related by 
L = XW. The remarkable feature of Little's Law is that 
it holds regardless of the specific forms of the arrival and 
departure processes. When applied to stochastic gene 
expression models, this implies that the processes lead- 
ing to mRNA/protein can be arbitrary, e.g. including 
gestation and senescence effects. 

We now apply Little's Law to derive an equation re- 
lating mRNA burst and steady-state distribution means. 
The arrival rate of mRNA bursts is driven by an arbitrary 
stochastic process with average arrival rate fc;,. The de- 
cay process of mRNA is also assumed to be driven by 
an arbitrary stochastic process with average decay time 
Tm- Employing Little's Law [l8|, we obtain a relation 
between the mean mRNA burst size nib and the average 
number of mRNAs in the steady state: 



(m) = At,„, 



(5) 



where A is average arrival rate of the mRNAs, which is 
given by 



A = rribkb 



(6) 



Hence, we derive that the steady-state distribution mean 
for mRNAs is related to the mean mRNA burst size by 



(to) = mbhTrr 



(7) 



Both the mean steady-state mRNA levels {{m)) and 
the mean mRNA lifetime (Tm) can be determined exper- 
imentally using standard procedures. Eq.[7] implies that 
the degree of transcriptional bursting can then be deter- 
mined by estimating the mean burst arrival rate (kb), 



which can be done using single-molecule approaches. 
Such a procedure was used in Ref. lj| to estimate the 



degree of transcriptional bursting, with the assumption 
of constant mRNA arrival rates and decay rates. Eq.[7]in- 
dicates that, even if this is not the case and arbitrary ges- 
tation and sensescence periods are considered, the above 
procedure remains a valid approach to determine the de- 
gree of transcriptional bursting mb- Alternatively, since 
the above relation is valid for arbitrary stochastic pro- 
cesses governing mRNA arrival and decay, it can serve as 
a useful consistency check for different experimental ap- 
proaches for quantifying mRNA burst and steady-state 
distributions. 

Using Little's Law we can also relate the steady-state 
protein distribution mean to the burst mean following 
similar logic. Since the average arrival rate of proteins is 
given by {m)kp, we derive 



-is — \ITT-}kpTp, 



(8) 



where kp and t~^ are average synthesis and decay rates 
of the proteins. The above equation can be recast in 
terms of the mean number of proteins produced in a sin- 
gle burst Pb (which is related to the mRNA burst distri- 
bution mean by Pb = nibkpTm)- Since the mean arrival 
rate of proteins is given by kbPb, we have 



^s — '^b^b'^'p 



(9) 



It is noteworthy that this simple relation is valid for ar- 
bitrary gestation and senescence waiting-time distribu- 
tions. It establishes that the mean steady-state protein 
level only depends on the average protein arrival and 
degradation rates and is independent of the higher mo- 
ments of the corresponding waiting-time distributions. 
Thus, it explains the observation in Ref [11[ that gesta- 
tion and senescence do not affect the average susceptibil- 
ity to changes in parameters. 

Another important consequence of Eq.|9] is that pro- 
cesses that alter the burst distribution mean without af- 
fecting protein degradation times or burst arrival times 
will produce a proportionate change in the steady-state 
distribution mean. Thus, regulatory interactions which 
are sensitive to the degree of transcriptional bursting and 
alter protein burst distributions will produce proportion- 
ate changes in protein steady-state distribution means. 
This, in turn, suggests the possibility of obtaining signa- 
tures of transcriptional bursting by observing changes in 
steady-state protein distribution means upon regulation. 
To explore this possibility, let us consider how regulation 
by small RNAs modulates protein burst distributions. 



Regulation by small RNAs 

We consider regulation by small RNAs (sRNAs) based 
on a coarse-grained model (Fig. [1} studied previously 
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FIG. 1. The kinetic scheme for regulation of protein produc- 
tion by small RNAs with coupled degradation rate 7. 



I9I42I I which apphes to sRNAs that regulate mRNA tar- 
gets stoichiometrically due to coupled degradation [22| . 
Synthesis of sRNAs is taken to be a Poisson process with 
constant rate kg and the sRNA degradation rate is also 
taken as constant (/is) in the following analysis. The pa- 
rameter 7 controls mutual degradation of mRNAs inter- 
acting with sRNAs. As in the previous section, mRNAs 
are created in bursts, with the average rate of arrival for 
bursts given by fcf,. If ki,^s << 1, i.e. if the sRNA life- 
time is small compared to the mean arrival time between 
bursts, the distribution of sRNAs prior to a mRNA burst 
can be approximated by the steady-state distribution of 
sRNAs in the absence of mRNAs. Given this approxima- 
tion, we wish to derive expressions for the protein burst 
distribution in the presence of sRNAs. This is, in general, 
analytically intractable. However by employing further 
approximations which are valid for a range of parameters 
we can obtain analytical expressions for the burst dis- 
tribution. Specifically, we assume that synthesis of new 
sRNAs during a burst can be ignored, i.e. no new sRNAs 
are created in the time interval between mRNA creation 
and decay. Furthermore, we consider jTm ^ 1 such that 
mRNA degradation in the presence of sRNAs is assumed 
to occur due to mutual degradation with a sRNA rather 
than natural decay with average rate fim = —■ Given 
that these approximations are valid, a simple analytic 
expression for the mean regulated protein levels can be 
obtained as a function of mean sRNA levels as shown 
below. 



As indicated in the reaction scheme in Fig HI a pair 
of molecules of mRNA and sRNA can combine and be 
degraded rapidly with rate 7. We first consider the limit 
7 ^' 00. In this case, regulation by sRNA results in an 
instantaneous modification of the distribution of mRNAs 
just after the burst. The mRNA burst distribution prior 
to interaction with sRNAs is given by irynim) (Eq.[2]). 
The modified mRNA burst distribution after interaction 



with sRNAs (7rm("T-)) is given by 

CO 

7r^{m) = ^p{n)7T^{m + n), m > 1, (10) 

n=0 

where p{n) is the probability of finding n sRNA molecules 
at the time of burst. Any burst of m mRNA molecules 
instantly becomes an effective burst of m — n mRNA 
molecules (for m > n) due to coupled degradation with 
n sRNAs. Urn < n, the mRNA burst after the regulation 
will be effectively an 'empty' burst. The probability of 
an empty burst is given by 



7rm(0) = 1 - X! '^m('^)- 



(11) 



771 — 1 



Since the unregulated mRNA burst distribution is ge- 
ometric (with parameter p, say), we derive 



Ti'niim) = 7rin(m)^(l - pY p{r 



71=0 



= G(l-p)^,„(m), m>l, (12) 

where G{1 — p) is the generating function of sRNA 
probability distribution p{n), evaluated at the point 1— p. 
Using Eq.[Tl]we derive 



^„,(0) = l-G(l-p), 



(13) 



The regulated mRNA burst distribution is thus a con- 
ditional geometric distribution as in the unregulated case, 
but with modified average arrival rate fc,„ ~ G{1 —p)km- 
This is because 1 — 7rin(0) — G(l — p) is the probabil- 
ity that the regulated burst results in atleast 1 mRNA. 
Therefore, the average number of mRNAs in the steady 
state for the regulated case is given by (according to the 
equation Eq.[7]) 



(to) = rribkri 



G(l-p)(m). 



(14) 



We denote by Ph^ris) the mean burst size for proteins 
in the presence of sRNAs, where Ug = kg/ lis- Using the 
equations Eq. (|8ll4p for mRNA's steady state average, we 
derive in the limit of fast coupled degradation (7 — > 00) 



Ps{ns)^G{l-p)Ps{0). 



(15) 



Taking the sRNA distribution prior to the burst {p{n)) 
to be a Poisson distribution with mean ris = ks/p^s, and 
given that the mean mRNA burst size is given by mi, = 
1/p (Eq.|4]), we derive 



Pbjris) 
Pb{0) 



a/nib 



(16) 



Thus, if the burst mean (P^) is determined along with n^, 
the above relation determines mt, and hence the degree 
of transcriptional bursting. Eq.[5] further implies that 
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FIG. 2. The relative error rj of the estimated mt value (actual 
value rrib = 10) as a function of mean sRNA levels n^. The 
estimate from the flat portions of the curves is within 15% 
error for the parameters shown with -^ = 20. The results 
shown are for the case of exponential waiting-time distribu- 
tions for mRNA arrival and decay, comparable relative errors 
were obtained for more general waiting-time distributions. 



distinguish between the Poisson and Telegraph scenarios. 
The parameter ranges for validity of the above analy- 
sis are accessible experimentally based on previous work, 
e.g. high values of 7 relative to the natural degrada- 
tion rate ^m are expected for the sRNA RhyB [2j], and 
mRNA burst arrival rates which are small compared to 
the mRNA degradation rate have also been reported |l4| . 
Finally, we note that it would be of interest to apply the 
preceding analysis to systems which show high degree of 
transcriptional bursting primarily arising from random 
activation and inactivation of the promoter state j2S | . In 
particular, it was observed 25| that increasing concentra- 
tions of a transcriptional activator resulted in increasing 
the mean burst size rather than affecting the burst fre- 
quency. Since the procedure proposed in this work is an 
independent approach to determine the burst mean, it 
would be of interest to further analyze the above system 
using the analysis proposed in the current work. 



the ratio of protein steady-state means for regulated to 
unregulated cases ( p (q\ ) is equal to the corresponding 
ratio for the burst means in Eg. 1161 This in turn im- 
plies that the mean transcriptional burst size mi, can be 
determined by considering changes in mean steady-state 
protein levels. Taken together, these results provide a 
novel procedure for determining mi,. 

The proposed procedure has been computationally val- 
idated for a range of parameters using stochastic sim- 
ulations (Fig. 2). Specifically, we set up simulations 
based on the standard Gillespie algorithm [23| wherein 
the waiting-time for the next reaction is drawn from 
an exponential distribution. To consider effects such 
as mRNA senescence, we model mRNA degradation as 
a multi-step process, wherein the waiting-time distribu- 
tion for each step is drawn from an exponential distribu- 
tion such that the degradation time for mRNAs follows a 
Gamma distribution (see Appendix). Similarly, mRNA 
arrival was simulated as a multi-step process with gamma 
waiting-time distribution between mRNA arrival bursts. 
The output from the simulations is the mean steady-state 
protein levels as a function of the mean sRNA levels (ris), 
where the mean sRNA levels are varied by increasing the 
sRNA creation rate kg. Provided that the system pa- 
rameters are consistent with the following constraints: 
7 S> fJ-m ^ fJ-s, simulations indicate that the transcrip- 
tional burst size rub can be predicted with reasonable 
accuracy from the ratio of measured protein steady-state 
means for regulated to unregulated cases as discussed 
above. The errors in the estimate for mt using the above 
procedure are related to the validity of the approxima- 
tions made and are discussed further in the Appendix. 
Provided that regulatory small RNAs can be designed 
with parameters subject to the constraints noted, the rel- 
ative error in estimating mb is small and thus we can de- 
termine the degree of transcriptional bursting and clearly 



SUMMARY AND CONCLUSIONS 

In summary, we have considered a generalized model 
of gene expression with bursty production of mRNAs 
and proteins. Since very different stochastic processes 
can lead to steady-state distributions that are experi- 
mentally indistinguishable, the degree of transcriptional 
bursting cannot be inferred from steady-state protein 
distributions. In light of this, it has been argued that 
determination of transcriptional bursting requires dy- 
namic measurements of mRNA molecules in single cells 
111 Il2l . |26[ . In this work, we have derived exact relations 
connecting mRNA/protein burst and steady-state distri- 
bution means which are valid for arbitrary gestation and 
senescence waiting-time distributions. We further ana- 
lyzed how protein burst distributions are modified due 
to regulation by small RNAs for a range of parameters. 
Our analysis computationally demonstrates an alterna- 
tive procedure for quantifying transcriptional bursting, 
which involves measurements of changes in mean pro- 
tein steady-state levels induced by interactions with small 
RNAs. The strategy presented can also be applied to a 
broader classes of biological networks whose analysis re- 
quires inference of internal variables from observations 
at higher levels. An alternative strategy to direct mea- 
surements of internal variables is to discriminate differ- 
ent possibilities for the internal variables by coupling to 
a controlled external interaction. 



APPENDIX 
Finite 7 corrections 

The analysis in the main text considered the limit 
(7 — >■ 00) and we now consider corrections due to finite 



7 values. Lets take a more detailed look at the protein 
production process during the burst. We denote the du- 
ration from the beginning of the burst to the time when 
sRNA or mRNA number first reaches zero as the first 
stage of the burst. If the mRNAs outnumber the sRNAs, 
excess mRNAs will be left after the coupled degradation 
and evolve accordingly. We call the duration from 
this point to the time when all mRNAs are degraded 
as stage two of the burst. In the case that 7 — !■ oo, 
the duration of stage one will be zero and all proteins 
are produced in stage two of the burst. However, for 
finite 7 value, one has to take into account proteins 
that have been synthesized during stage one of the burst. 




100 



In order to estimate the amount of the proteins 
produced on average from mRNAs that are degraded by 
sRNAs (stage one), we observe first that the minimal 
degradation rate of a single mRNA in this process is 7. 
This is because at least one sRNA should be present to 
ensure coupled degradation. Second, the total amount 
of mRNAs in the originating burst is greater than or 
equal to number of mRNAs degraded by sRNAs (since 
some mRNAs may decay naturally). 

Hence, we can employ formula Eq.[7] in order to es- 
timate contribution of the mRNAs decaying in coupled 
degradation process to overall steady state level. The 
upper bound of this contribution is given by 
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(17) 



Here we replaced the rate r„^ in the Eq.[7] by the 
minimal rate 7 in order to estimate the upper limit. 

Now we can use the expression Eq.|8] to get the up- 
per bound of the proteins produced on average from the 
mRNAs during the coupled degradation process 



SPsirig) ~ SmkpTp ^ n. 



7 



f^p'^p- 



Hence, the overall ratio of regulated to unregulated 
mean steady state levels of proteins is bounded as 
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(19) 



which is independent of protein's synthesis rate kp. 
Therefore, if coupled degradation process is much faster 
than natural mRNA decay, Tmj ^ 1, we obtain 6R — >■ 
which is validated by simulations. As we can see in 
FigQ, when t^J > 10 the proteins produced during 
stage one of the burst can be neglected and the result is 
almost the same as when 7 — > 00. Finally we note that 
recent studies [2J| have shown that a well-studied bacte- 
rial small RNA (RhyB) does induce rapid degradation of 
target mRNAs consistent with the condition Tmj 3> 1. 



FIG. 3. The relative error 77 of the estimated mt derived from 
changes in mean protein burst levels. For different Us and mi, 
values, the error is negligible when 7 > lOfim- 



Waiting-time distribution for multi-step processes 



Previous work [ll| on gestation and senescence effects 
in mRNA/protein production and decay considered ex- 
tensions of the single-step Poisson process to multi-step 
processes. For the simplest case, the corresponding wait- 
ing time distribution is a Gamma distribution as derived 
below. Consider a multi-step process, consisting of n 
steps such that each step is completed with rate k. Let T 
denote the random variable corresponding to the waiting- 
time for the process to finish and let Ti be the random 
variable corresponding to the waiting-time for the i 
step. Thus we have T — J^i'^i^ i-^- ^ i^ ^he sum of 
n identical independent random variables. Correspond- 
ingly the Laplace transform of the probability distribu- 
tion for T (denoted by F{s) say) is given by the product 
of n Laplace transforms of the exponential distribution. 
(18) The exponential waiting-time distribution for the i^^ step 
is given by ke~^* with corresponding Laplace transform 
-^r^- Correspondingly we have F{s) = (j^r^)", and in- 
verting the Laplace transform we obtain that the waiting- 
time distribution for the multi-step process is given by 
the Gamma distribution: k{kt)'"'^^ (k-iv. • 
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